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METHOD OF DERIVING ORBITAL PERTURBING PARAMETERS 


FROM ONBOARD OPTICAL MEASUREMENTS OF 
AN EJECTED PROBE OR A NATURAL SATELLITE 

INTRODUCTION 

The principal objective of this analytical investigation is to 
evaluate the usefulness of a number of sightings taken from an orbiting 
spacecraft to another (secondary) satellite. The sightings determine 
the direction of the secondary satellite relative to the fixed stars 
by angle measurements only. The sightings are assumed to be made 
intermittently over periods ranging from a fraction of one orbit to 
many days. The secondary satellite may be a probe which is ejected 
from the spacecraft itself or a natural satellite. 

This program has been exploratory in nature and, in general, is 
concerned with the possibility of deriving a knowledge of forces which 
perturb the orbit. From a theoretical viewpoint the feasibility of 
navigational systems of this type is already established. What is not 
known, however, is the effect of measurement errors on the accuracy of 
various perturbing parameters. Further, the determination of the 
perturbing parameters and the solution of the navigation problem con- 
stitutes a single problem for which solutions must be obtained simul- 
taneously. 

The study is conducted primarily with respect to spacecraft orbits 
about Mars, Venus, the Moon, and Earth. Some results are shown for 
Saturn. The study of Mars utilizes various combinations of one or 
both of its natural satellites and an ejected probe as secondaries. 
However, for orbits about Venus, the Moon, and Earth, the secondary is 
a single ejected probe. The speed of ejection is very small as compared 
with orbital speeds. 
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PROGRAM OBJECTIVES 


Figure 1 indicates the geometry of the sightings and the parameters 
studied for the case of orbits about Mars and Venus. These parameters arise 
from gravitational forces. In a following study concerning the moon and the 
earth, a more complicated potential will be assumed for the moon, and non- 
gravitational forces will be considered for the earth. 


As a starting point, effort has been concentrated on the problem of deter- 
mining the higher order terms in the gravitational potential of Mars. Since 
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Figure 1: Perturbing Parameter Study for Mars and Venus 
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Mars has two readily viewable natural satellites (Phobos and Deimos) , it was 
first decided to base the investigation on angular observations of these bodies 
rather than to assume that an artificial probe was ejected. This, in turn, 
necessitated a prior navigational study involving 12 and 18 unknowns to deter- 
mine the basic characteristics of such systems from the point of view of error 
propagation. As a result, it was found that systems in which our own satellite 
orbit was unknown as well as the orbits of Phobos and/or Deimos were indeed 
promising from a navigational point of view. Although the 18 unknown problem 
had satisfactory properties, it was decided that observations of a single 
orbiting body would be operationally simpler and thus the 12 unknown problem 
(our own satellite and Phobos) would form the starting point for Mars perturbing 
parameter studies. 

After completing the Mars perturbing parameter study based on the use of 
Phobos as a secondary, consideration was given to including sightings on an 
ejected probe. This was done and proved to yield more accurate results, as will 
be shown. 

The next phase of the investigation involved the determination of the 
gravitational potential of Venus. Since Venus has no known natural satellites, 
an ejected test probe of assumed ejection speed but unknown direction was used. 
The effect of an error in the assumed speed is also studied. The degree of 
interest in the various perturbing parameters is different for Venus than for 
Mars. Because Venus is rotating very slowly (rotational period of 247 days 
retrograde), it is probably the most nearly spherical of the planets. It turns 
out that considerable difficulty is encountered if the direction of the axis 
of dynamic asymmetry is treated as an unknown in the total problem. 

In the moon study three separate orbits were chosen, each of which had 
an eccentricity of 0.1; the secondary being a probe ejected at 4 km/hr. In 
this case, problem emphasis is quite different from that of orbits about Mars 
or Venus for the following reasons: 

(1) The moon has a complicated potential in that tesseral harmonics 
must be considered. 
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(2) Earth-based radar measurements are more effective in the case of 
the moon or the earth than Mars or Venus. Hence, our method must 
be compared with stronger competition. 

As many as 10 unknown harmonics coefficients were chosen. In general, our error 
analysis yielded accuracies comparable to those obtained from Lunar Orbiter data. 

Finally, low and high altitude Earth orbits were investigated. For this 
case the potential was assumed known, and unknown drag and radiation pressure 
terms were allowed. It was found that if a small error in the speed of probe 
ejection were made, a significant error in these parameters would be incurred. 

In order to overcome such errors, a geometry was studied in which the probe 
was simply released; drag and radiation pressure were then allowed to separate 
the probe and the spacecraft. 

It is one of the principal objects of the investigation to determine the 
number of orbits over which measurements must be made in order to determine the 
values of the parameters for various degrees of accuracy. The criterion which 
is used in determining what constitutes an acceptable degree of accuracy is the 
accuracy with which the parameter in question is known at present. Thus, in 
the case of the earth, accuracies must be very great, whereas in the case of 
Mars and Venus, where large uncertainties exist, acceptable accuracies need not 
be as great. 


Measurement Techniques 

In performing the present study emphasis has been placed on the analytical 
problem and questions of measurement are subsidiary. However, it is desirable 
that mention be made of possible methods of measuring the position of the target 
relative to the celestial coordinate system as well as the problem of detecting 
the target under various viewing conditions. 

It is assumed that the angular coordinates of the target are measured to an 
accuracy of ten seconds of arc for each sighting. For most of the problems 
simulated it was assumed that about 13 sightings were made per orbit. If the 
view of the target was blocked by the central body, then that sighting was 
omitted. Otherwise, the sightings were uniformly spaced. 
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In all cases the basic coordinate frame is that which is supplied by the 
catalogued position of the stars. Several types of detection instruments could 
be used to measure the direction of the target in this frame, e.g. a star 
tracker or an optical scanner. A review article describing optical scanners 
and their applications is given by reference M- 

Planetary Constants 

The specific technique used in the study of the perturbing parameters 
depends on the physical constants of the planetary body in question. Tables 
I, II, and III contain relevant data. The quantities in parentheses behind 
each number indicate the source. 

Table I gives some physical constants for Mars, Venus, Moon, and the 
Earth. As can be seen from this table, two of the bodies, the Earth and Mars, 
rotate sufficiently fast to have a measurable flattening which, in turn, permits 
the computation of the right ascension and declination of their spin axes as 
well as several terms in the gravitational potential. (As is shown in Table II, 
this is true even to a greater degree for Jupiter and Saturn.) The other two 
bodies, the Moon and Venus, are sufficiently close to spherical in shape that 
a more general representation of their irregular shape is necessary, and 
further, the gravitational perturbations only weakly define the orientation 
of the spin axis. 

In studying Mars, we utilize the natural satellites, Phobos and/or Deimos . 
Some of the characteristics of their orbits are shown in Table III. 
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Table Is Physical Constants for Mars, Venus, Moon, and Earth 


Constant 

— 

Mars 

Venus 

Moon 

Earth 

Ratio of body mass 
to Earth mass 

0.10730 (7) 

.81500362 + .123 x 10" 5 
(2) 

1/(81.3004) + .0007) 
(2) 

1.0 

Radii (km) 

R = 3388 + 26 
e — 

18 < R - R <36 

R = R = 6100 + 50 
e p - 

R x = 1738.549 + .061 
R 0 = 1738.189 + .061 

R = 6378.15 + .05 
e - 

R = 6356.63 + .05 


6 P (7) 

(6) 

R, - 1737.470 + .061 
(9) 

P 

Reciprocal 

Flattening 

150 + 50 (7) 

Large (nearly spherical) 


298.30 + .05 (3) 

Potential 

J 3 

J 4 

J 5 

J 6 

2.011 x 10" 3 (11) 


2.07 x 10" 4 

- 0.4461 x 10” 4 

- 0.2089 x 10“ 4 (4) 

(1082.28 + .3) x 10“ 6 

(- 2.3 + .2) x 10“ 6 

(- 2.12 + .5) x 10“ 6 

(- .2 + .1) x 10“ 6 

(1 + .8) x 10‘ 6 (3) 

or (9) 

Rotational Period 

24 h 37 m 22.6689 S 
(7) 

274 + 4 days retro- 
grade (6) 

27.32166 days (9) 

23 h 56 m 4 s 

Spin Axis r.a. 

Direction 

316.55° + .00675° 
(t - ^905) 

272.75° + 4° 

A A 

~ - 91° 

undefined 

dec. 

52.85° + .00346° 
(t - 1905) (7) 

71.50 + 4 (6) 

~ 88.5° * 

90° 


* On Q page 51 of the American Ephemeris and Nautical Almanac i and Q 1 are tabulated. Declination = 
90° - i, right ascension » 90° -"h 1 . 




Table II: Physical Constants for Jupiter and Saturn, From [3] 


Constant 

Jupiter 

Saturn 

Ratio of Sun's mass 

1047.355 + .065 

3,500 + 3 

to body mass 



Reciprocal Flattening 

15.2 + .1 

10.2 

Major Radii (km) R 1 

71,375 + 50 

60,500 + 50 

*2 

66,680 + 50 

54,550 + 50 

Potential J 2 

.0148 + .0001 

.0161 + .00001 

J 4 

- .00035 + .00005 

- .000865 + .000026 


Table III; Constants of Martian Satellites 


Constant 

Phobos 


Deimos 


Mass 

Diameter (miles) 

10 

(7) 

5 

(7) 

Characteristics of Orbit 





Orbital Period (days) 

0.318910 

(5) 

1.262441 

(5) 

Semi -axis (sec at 1 AU) 

12.895 

(7) 

32.389 

(7) 

Ellipticity 

0.0210 

(7) 

0.0028 

(7) 

Orbital plane 

see ref. 

(7) 

see ref. 

(7) 

Regression of ascending 

158.5 

(7) 

6.5 

(7) 

node (deg./yr.) 

** 


■3 


Mean distance from 

9.4 x 10 

(5) 

23.5 x 10 

(S'* 

planet (km) 





LIST OF MAJOR SYMBOLS 


The following is a list of major symbols used in the sections 
ANALYTICAL STATEMENT OF THE PROBLEM and METHOD OF SOLUTION, 
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Rli(t) 

*2i00 


celestial coordinate system 


inertial coordinate system with parallel to planet's axis of 
symmetry 

th 

coordinate system which is associated with the i Kepler arc which 
approximates the spacecraft’s orbit for (i - 1) T < t < ir 

th 

coordinate system which is associated with the i Kepler arc which 
approximates the secondary's orbit for (i - 1) t < t < iT 


t = time 


t = fixed time length 

- angles which define orientation of with respect to 
R^(t) = position vector of spacecraft 
(t) = position vector of secondary 

k = vector of unknown parameters (constants) 


XjOO = 


( R j (t A 
V (t) / 


j = 1, 2 


e = declination of secondary's position as viewed from spacecraft 
a = right ascension 

*■’" th 

Rli K (t) = approximate position vector of spacecraft along i Kepler arc 

th 

RzikOO “ approximate position vector of secondary along i Kepler arc 

6X (t) (t) - X. ft), (i - 1) T < t < iT 

3 J 

e , a = angles which define direction of probe ejection with respect to 

r j 5 S j ,Z j = com P onerits expressed in R^(t) 

£v= ejection speed 

In general. 


U* = transpose of U 

. = fL 

dt 

a = implies u • u = 1 


6u = error in u 

a(u) = standard deviation of a random variable u 
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ANALYTICAL STATEMENT OF THE PROBLEM 


In general, if a second satellite is sighted from a primary satellite, 
the time-path observed is a function of: 

(i) the position and velocity of the primary at t = 0, 

(ii) the position and velocity of the secondary at t = 0. 

(iii) the force field in which the two satellites are traveling. 

Since the observed time-path is dependent upon these three factors, it 
may be possible to indirectly measure these factors by making measurements at 
known discrete times of the secondary* s position with respect to the primary’s. 

Equations of Motion 


Since the motions of the primary and secondary satellites are governed by 
a sixth order system c£ differential equations, it is possible to write 


R x (t) = G(R 1 (t), R 1 (t), t, k) 
R^t) = G(R 2 (t), R 2 (t), t, k) 


( 1 ) 

( 2 ) 


where R^(t) - position of the primary as a function of time, 

R^Ct) = position of the probe as a function of time, and 
k = a vector of unknown parameters (constants). 

Instead of writing the system of equations as three second-order equations, it 
is more convenient to consider six first-order equations. To this end, let 


X.(t) 


R.(t> 


a 6 x 1 vector , i = 1, 2 . 


Thus, Equations (1) and (2) may be written 


X, 



^ > k) 9 i 1,2. 


We will assume that F may be written as 
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_ n _ _ 

F(X, t, k) = F (X) + 2 k P CX, t) 
K j=l 3 3 


where 
each Pj 


is the vector which arises from the inverse square force field, and 
is a known function of its argument. Hence, (1) and (2) take the form 


X. = 


w 


n 

+ 2 

j=l 


k. P.(X., 

J J i 


t), i = 1, 2. 


(3) 


Constraint Equation 


If measurements of the position of the secondary with respect to the 
primary are taken, these measurements imply a constraint on the motion of the 
two satellites. In this report, we will assume the measured quantities are 
the right ascension and declination of the vector from the primary to the 
secondary at known discrete times. That is. 


R 2 (t p ) - R^tp) = (cos e cos a 


i^ + cos e sin a 


j ^ + sin e 


(4) 


where 



m. = ivy - vyi 

e = declination of (measured) 

a = right ascension of R^ - R^ (measured) 

k = unit vectors associated with the celestial coordinate system, i 
in the celestial equatorial plane and in the direction of the 
first point of Aries, k^ in the direction of the North Pole, j 
completes the right-handed coordinate system. 


Equation (4) requires that R^ - R^ be written with components in the celes- 
tial coordinate system (we will denote this system as S^). However, this co- 
ordinate system may not be convenient when considering the equations of motion 
(3). In order to minimize n in (3), a new coordinate system is chosen with 
one axis parallel to the planet's principal axis of largest moment of inertia. Two 
angles, and are required to define the orientation of this new. system, 

S^, with respect to the celestial system S^. Hence, 

i^ rotation + 90° about 

— k^ rotation about f 2 = i^. 
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Hence, (4) becomes 


d( tp ) = ^Ax(t p ) 


where 



x = R 2 (t p ) - Rl (t p ) 

= Xi £3 + x 2 5 3 + x 3 k 3 


IT 


2 2 
+ x 2 + x 3 


( 5 ) 


ments 


Equation (5) then is the basic constraint equation implied by the measure- 
. This equation yields two independent relations which constrain x(t ). 


The Problem 


The problem may now be stated as follows: Find initial conditions 

X^(0), X 2 (0), and force field parameters k ^ , j = 1, 2, . .., n, which imply 

X^(t) and X 2 (t) as solutions to (3) such that X^(t) and x 2 (t) satisfy (5) at 

the discrete times t • 

P 

In some cases, the angles ^ and must also be treated as unknowns. 

This is particularly true for the case of orbits about Venus. Also, if the 
secondary satellite is an ejected probe, then at the time of ejection. 

A complication to the analytical problem is that the number of measure- 
ments is expected to be much greater than the number of unknowns. Hence, 
because of various errors, no X^(0), X 2 (0), and k will exist such that all 
constraint equations (5) can be satisfied. Thus, some "best" value of the un- 
knowns must be found. 
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METHOD OF SOLUTION 


There are many methods of attacking the formulated problem, but we will 
discuss only one particular approacji. The basic technique makes use of a 
linearization of the equations of motion and the constraint equation. Hie 
linear equations are then solved to yield improved values of the unknowns. An 
iterative computational procedure and an initial guess are required. Linearization 
of the constraint equation is as follows: 


The total differential of (5) yields 

6d = — A[6x - - + — 6A x ; 

M. ' 

but, ^fip, - x • 6x, since V? = x • x. 

So, p,fid = A[fix - — — x] + 6A x 
M- 

where 


6d = 


cos e sin a 6 a - sin e cos a fie 

cos e cos a 6 a - sin e sin a fie 

cos e fie 


So, p,6a cos e = (- sin a ^ + cos a A 2 )[6x - X x ] + (“ sin a 6A 1 + cos a 

^ • 6A 2 ) x 


where 


A. = k tk row of A, 
k th 

fiA^ = k row of fiA. 


Also, 


M-fie cos e = A^Cfix - - x ] + 6A^ x. 


But, 
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k tli 

where A = k column of A. 
Hence , 



Here, 


l±j_ “ - sinCS! - a) - cos^ - a) cos § 2 x 2 + cos^ - a) sin § 2 x 3 , 
M- 2 = " a) (sin § 2 x 2 + cos § 2 x 3 ), 

= sec e (cos § 2 x 2 - sin § 2 x 3 ). 


Equation (6) is two independent relations in §a, 6e, 6g , 6§ 2> and 
6x(t p ) = 6R 2 (t p ) - SR^tp). 
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Linearization of the Equations of Motion 


For the present, let us drop the subscript i from Equation (3). Now let 

X(t) = Xj^t) + 6X(t) (7) 

k = ko + 6k 


where X^(t) is a reference trajectory (known) of possibly unconnected Kepler 
arcs . Thus , 


' W 


except for a finite set of points. Substituting (7) into (3) we obtain 


6X = 6X + .. 

„ _ n 

= VF k (X k ) 6X + 2 

j=l 


n __ _ _ _ 

+ k j ( P j (X K» fc > « x + ...) 

(ko j + 8kj) t). 


( 8 ) 


Thus, to a first approximation, 6X satisfies a linear non-homogeneous differen 
tial equation. We may now write the solution to (8) as (Coddington and Levin- 
son [8j) 


t 

6X(t) = §(t, t 0 )[6X(t 0 ) + J 

to 


$ _1 (s. 


n 

t) 2 

J-l 


k^ P j (X K (s), s)dsj 


(9) 


where 

I = VF k (X k ) $ and $(t Q , t Q ) - I. 

The matrix § in (9) is commonly called the transition matrix in the term- 
inology of orbital theory. It is the fundamental solution matrix of (8) in 
the terminology of differential equation theory. 

An explicit solution of $ is given by Kochi [LO] . To utilize his solution 
the components of 6X(t 0 ) and 6X(t) must be expressed in preferred coordinate 
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systems. This point will be considered next. In any event, (9) gives a 
simple method of calculating an approximation to the deviation from a single 
Kepler orbit given 6X(t 0 ) and k. 

In deriving an explicit expression for $, Kochi introduced a coordinate 
we will call R(t) which moves with the Kepler reference position R^(t). bet 
R(t) have associated unit vectors r, s, z such that r = R^(t ) | , s is 

in the plane of the Kepler orbit so that R^(t) • s > 0, and z completes the 
right-handed system as shown in Figure 3. 


APPROXIMATE POSITION 



Figure 3: The Moving Coordinate System R(t) 

With this definition, we may write 
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where 


B(t) 


CQ CtU 
SQ CUJ 


- sQ ci scu 1 
+ cQ ci Scjo 1 
si su ) 1 


- cCl su) f - sQ ci ecu' 

- sQ sod 1 + cQ ci ecu' 

si ecu 1 



where c0 = cos 0 
s0 = sin 9 

i = inclination of Kepler orbit 
Q = longitude of line of nodes 
u; = argument of perigee 
cu' = w + v (t) 

where 


v(t) = true anomaly. 


Since we will use $ as developed by Kochi, in (9) 6X(t Q ) must have posi- 
tion and velocity components expressed in R(t 0 ) while 6X(t) must have these 
components expressed in R(t). 


Since the perturbing forces are to computed along the reference 


orbit, no single Kepler orbit can be expected to yield accurate results. We 
shall choose discontinuous arcs of Kepler orbits as a reference orbit as shown 
in Figure 4. Each arc is of time duration, r. 


In general, we have for (i - 1) t < t < It 

6X(t) = §.(t, (i - 1) t) [6X«i - 1)T) + 2^ (k oj + akj) F.-jCt)] 

where x iK ( fc ) is the position and velocity along the i Kepler arc, 

$. is evaluated using this arc as a reference. 


? ij (t) = J (S ’ (i " 1)T) P j (X iK (s)> s) 


ds 


(i-l)T 


6X((i - 1)t) must have components in R^((i - 1)t) 
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ARC OF REFERENCE ORBIT 



Figure 4: The True Orbit, Reference Orbit, and Computed Orbit 


P.(X.^(s), s) must have components in R. (s) 
j IK 1 

6X(t) must have components in R^(t). 

But we will choose to define 6X(1t) so that. 

6X(iT) = i|r.[~6X((i - 1)t) +2 6k. E. .1 

1L j = 1 J 1JJ 


where = §^(iT, (i - 1)T) 


E. . = F. . (iT). 
ij xj 


Thus, to generate the i + 1 Kepler arc, we use the conditions 
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X i+l,K (lT) X iK (iT ^ + j f 1 ^ j 



( 10 ) 


In (10) X and X. 1 v must be expressed with components in S~. After com- 
— i+l, K. 

puting X. (iT), we may find the orbital elements of the i+1 Kepler arc, 

l+l, K 

and hence $_£ + ^(t, i T )* So, for iT < t < (i + 1)t , 

6X(t) = $ i+1 (t, iT) [6X(ix) + 2^ k j F i+1>j (t)] 

B| +1 (iT) B.(iT) 0 

0 B i+l (iT) B i (iT) y 

x [ 6 x (i - i) T ) + ^ ekjl..] +£ y^ct) 


the matrix B| + ^(iT) B^(iT) differs from the 3x3 identity matrix, but its 
effect is second order. Hence, for iT < t < (i + 1 )t 


6X(t) = * i+1 (t, i T ) 


6X(t) = $. +1 (t, 1 t)[m. 6X(0) + ^ 6k.. 


.(V. . + F. _ .(t)) 
’ xj x+1, J " 




iT < t < (i + 1 )t 


6X((i + 1 )t) = M. +1 6 X (0) + ^ 6kj V. +1J 


(U) 


where 


■ *i E ij + *i *1-1 E i.i,j + ••• + *i *1-1 


tj Ey 


M 1 ' *1 * 1-1 


‘I' 


So that. 
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V i+l,j ^i+l^i+l.j + V ij^ 


M. = i|r M. . 
l+l l+l 1 


Equation (11) thus gives us the deviation in position and velocity of the 
reference position and velocity from the time quantities at t, iT < t < (i+l)T, 
in terms of these deviations at t = 0, ^ , and file. In this equation $X(t) has 
components in R i+ ^(t). 


Linear Equation in Unknowns 

Equation (6) is two independent equations expressing a relationship be- 
tween errors in the measured direction to the secondary satellite, errors in 
the two angles which define the planetary coordinate system with respect to 
the celestial coordinate system, and errors in the position of the secondary 
with respect to the primary. This equation is useless in itself, for if 
errors in the measured quantities are specified (6a and §e) the equations con- 
tain more unknowns than independent relations. The dynamics of the problem, 

however, do not allow the vectors 6R- to move at random. In fact, if 6R.(0), 
— i 

6R i (0), and k are specified, then 6R^(t) and 6R^(t) are determined for all 
time, t. This dependence is analytically stated by Equation (11). In this 
section, we will utilize Equations (6) and (11) to obtain a relationship be- 
tween 6a, 6e , 6§ 1# 6X^(0) and 6X 2 (0). 

A trivial transformation yields 

6R.(t) = K 6X.(t) 

where K = (I, 0), I being the 3x3 identity; 0 being the 3x3 zero. So, (6) 
becomes 

cos e \ _ _ 

J = E K(6X 2 (t p ) - 6X x (t p )) + F 
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In the above equation 6X^ and 6X^ must have components in S^. In order to 
utilize this equation, it is necessary to rewrite it so that SX^Ct^) has com- 
ponents in R^i(t^) while 6X^(t^) has components in R^(tp), where 


(i - 1) T < t < iT« 

- p 


Hence , 


6a cos e 


6e 


= E(B 2 .(t ) K 6X 2 (t ) - B u (t ) K 6X x (t )) + F 



So, Equation (11) yields 


" °2 + D 1 + 


6a cos e 
6 e 


| ( ^2’ " ^1’ ^21 ' ^11’ ^22 " ^12”“^2n " ^ln ,F) 


x 6Y (0) 


( 12 ) 


where 


(t> lT> £ ‘•J V,) ( V 

D 1 - E B li k ‘h <*• 

t 2 - « >> 21 «9 K * 21 (t p , (1 - l)r) 

lt 1 - £ B ll(tp ) K <i - Dr) M 2j . . 2 

* 2J ‘ E B 2i<V K *21 V (1 • ° T) <»2. i - 1, j + h, i, jV> 


v u ' E B u (t p )K *u <V (1 - « T > <»1. i - 1, J + *1, i, J<V> 


, a (14 + n) x 1 vector. 
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Thus, (12) is two equations in 6Y(0). If a sufficiently large number of 
measurements are specified we may solve for 6Y(0). After finding this vector 
6 X 2 (t) and 6 X^(t) may be found from (11)* Since (12) is only an approximate 
equation, an iterative computational process is necessary. 

Let us now suppose the secondary satellite is a probe which was ejected 
from the primary with a known speed but poorly known direction at t = 0 . Ihen, 

RjCO) = ^(0) , 


and 


r 2 co> = 


R^O) 


+ Av u 


where a v = speed of ejection (assumed known) and u defines the ejection direc- 
tion. We may choose the reference trajectories so that 


R 2 i K (°) r 11k (°) 


R 21 X^ = R 11 K^ + Av u 0 , where u 0 is an estimate of u. 


So, 6X 2 (0) = 6X 1 (0) + Av 



Now let u = cos e cos a 


i 3 


+ cos e sin a + sin. e 


Thus , 


6X 2 (0) = 6X^0) + 



<13) 


where 
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Q * Av 




and the components of 6 X 2 ( 0 ) and 6X^(0) are in S^. 


In (12), 6 X 2 (0) has components in 6X^(0) has components in 

Rll(0); (13) cannot be substituted into (12) without rewriting (12) such that 
6 X 2 ( 0 ) and 6X^(0) have components in the same coordinate system. To this end, 
let us rewrite (12) so that 6X 2 (0) and §X^(0) have components in S^. Hence, 


D 2 +D l + ! 


6 a cos e 
6 e 


.... \ n - TT ln , n«w)o*) 


where 


luT = 1 1 

w 2 2 


= 11 

w i 1 






B h<°> 


(2 x 6 ) matrix 


(2 x 6 ) matrix 


and 6X 2 (0) and 6X^0) have components in S^. Hence, (13) allows us to rewrite 
(14) as 



Ax 1 (0)\ 

$2 ~^l> ®21 " Tn* •**’ ^2n " ^ln’ F » P 2 Q ^ | 6k 

6§1 

l ^ I 

where ±s a 2 x 3 matrix defined by the last three columns ofllJ^. 


Other cases may be considered. For example, Av unknown; Av unknown; the 
time of ejection unknown; etc. 
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Summary of Equations 


We now list the equations as developed in the previous sections. 


1. We are given 

(a) x ^(°) and X 2 (0) in S 3 

(b) k Ql> k Q2 , k Qn 


(c) P^, P 2 , P^ as a function of spherical coordinates and with compo- 

nents in R(t) (for both trajectories) 

(d) times of observation t^, k = 1, 2, f 

( e ) ?2 
(f) t 

( 8 ) 



- cos cos § 2 

- sin cos § 2 

sin g 2 


cos sin § 2 

sin sin § 2 

cos ?2 


2. M 1q = M 20 = I (6 x 6), V 1Qj = V 2Q . =0 (6x1) 


3. X((i - 1)t) = X iK ((i “ !) T ) ( for both trajectories in S,,) 


4. Find orbital elements 


5. If t > iT 

k — 

(a) ^ = $ i (iT, (i - 1 )t) 

(b) M. = * i M 1 _ 1 

iT 

ijCiT) = f 4' 1 (s, (i - l)T)P i j(s) 

(i-l)T 

(d) V. . = t- (E. . + V. . .) 

ij *i v ij l-l, j' 

(e) C. = E k,, . lk. E. . 

i j=1 Y i ij 


(c) E. . = F. . 
ij ij 


ds 



(X. K (i(T)) + C.) 


B^iT) 


in S. 
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6 


. If < iT , ve obtain two linear equations in the unknowns. 
Use i Kepler arc to compute 

(a) mean anomaly at t^ 

(b) true anomaly 

Ml -<6 P, 

(c) r = — ^ = J 

j 1 + e. cos v # 1 + e. cos v. 

J J J J 


(d) 


(e) 


(f) 

(S) 

(h) 



x(t k ) = Xl (t k ) t 3 + x 2 (t k ) j 3 + x 3 (t k ) « 3 


/ cos e cos a 
I cos e sin a 
V sin a 


= — A x, determines sin e, cos e, sin a, cos a 



= - s in (g x - a) X x - cos^ - a) cos x 2 + cos^ - a) sin § 2 

P'2 = sin (§i - a ) sin 5 2 x 2 + sin ^l " a ) cos ?2 X 3 
cos § 2 sin | 2 

^4 cos e X 2 cos e X 3 
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(i) F_(t: k ) = J i:\s, (i - 1 )t) P^fs) ds 

(i-l)T 

3 = % 9 • • • 9 n 

11 = E B^t^K § 1 (t k , (i - l)r) M i l 

Vj - E W K »i<V (i ' 1 )T) ^i-l,j +? i, j (t k» 


(j) A k 
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F) 



2k x (14 + n) 


8. After all measurements the covariance matrix of the output errors is given by 


E(6Y T (0) 6Y(0)) = ^(aL)" 1 . 

PERTURBING FORCES 


In this study, the only perturbing forces acting on orbits about Mars, 
Venus, and the moon for which numerical results will be shown are those which 
arise from higher order terms in a gravitation potential. Let us now give a 
short discussion of the gravitation potential. 

Let (X, Y, Z) be an arbitrary set of right-handed orthogonal axis, and let 
(r, 0) be the spherical coordinates of a point, P, external to a body whose 

density is sectionally continuous (Figure 5). 


The potential of a unit mass at P due to the body may be written as 


V (r , *5, 



B 


D (P') dv 
|P' - P| 


(15) 
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X 


Figure 5: The Body and External Point, P 


where B = 
P' = 
P = 
G - 
D(P') = 
dv = 
Now, -W= 


region occupied by the body 
generic position within the body 
position of external point, P 
universal gravitational constant 
density of body at P* 
element of volume of P' , 

the gravitational force on a unit mass at P due to the body. 


It would be most convenient if the integrand of (15) could be written as 
f(P’) g( p ); for if this were the case, then, 

V(r, rf, 0) = -G g(P) JJJ f(P') dv 

B 

= G x g(P), 

where is a constant which depends only on the shape and mass distribution of 
the body. Hence, the force on a unit mass would be known except for a single 
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multiplicative constant. Unfortunately, the integrand of (15) cannot be written 
so simply, instead we must write it in the form 

CO 

£ f (P 1 ) g (P) . 
i=l 

So, 


V (r , <5, 0) = Z G. g (P) (16) 

i=l 


where each is independent of P and dependent only on the shape and mass 
distribution of the body. The exact form of (16) may be found in [9], and can 
be written as 


V(r, rf, e> = - 1 + 2 2 (■=*) P™(cos »5)(C 

L n=l m=0 \ / 


cos m 0 + S sin m 0 
nm nm 

(17) 


>] 


where C and S are constants which are dependent on the mass distribution 
nm nm 

of the body. 


m = total mass of body 

r = equatorial radius of body 
e 

G = universal gravitationa 1 constant 

P^(x) = associated Legendre function of the first kind of degree k and 
order j . 

A tremendous simplification of (17) occurs if the body is symmetric about 
the Z axis. If this is the case, then 


C = 0, m ^ 0 
nm 

S =0 


nm 

and (17) is customarily written in the form 


V(r , rf) 



where J. are constants (J, = -C ) 

k k ko 

P^(x) = Legendre polynomial of degree k 


p£(x) 
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Z) is the center of mass of the 

*>] ( 18 ) 

We will use this form in the study of Mars and Venus in this report. Recall 
that the validity of this form rests on two assumptions: 

(1) the Z axis from which the angle 4 > is measured is an axis of symmetry 
of the body , and 

(2) the distance r is measured from the center of mass of the body. 

The assumption that Mars and Venus have an axis of symmetry appears to be 
a good approximation, for there are some reasons to believe that these planets 
are more homogeneous than Earth. However, the direction of the axis of symmetry 
with respect to the celestial coordinate is not well known beforehand. Hence, 
if we choose to use form (18) we are forced to introduce two angles which de- 
fine the direction of the axis of symmetry with respect to the celestial system 
as unknowns. 

The assumption that the coordinate system may be located at the center of 
mass of the body introduces no significant difficulties and will be used here. 

Let us now investigate the change in the coefficients of (18) for the case 
in which the coordinate system chosen lies close to the preferred system. 

As before, let be the celestial coordinate system. Also, let be a 
system with associated unit vectors, i^, and such that *- s parallel 
to the axis of symmetry. Let and the azimuth and coelevation of 

with respect to (Figure 5). Then, the potential has the form (18). 

Now suppose we do not utilize the preferred system, S^, but instead write 
the potential in a system defined with respect to by angles 5^ + 6?^ and 
^2 + 6 ^2* Ifc may be s * lown that for sma H and 6^ 


In addition, if the origin of the axes (X, Y, 
body, the = 0, thus (17)becomes 


V (r , « - - ^ 1 - S j 
k=2 


W 


P k (cos 
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V (r , 0) 



where 

°kl • - J k 6 ?2^ 
S kl ' ' J k si ” ?2 


i> = coelevation of P in new system 
*0 = azimuth of P in new system. 


Thus, to a first approximation, Legendre functions of the first order must 
be included in the potential if a coordinate system slightly removed from a 
preferred system is utilized to describe the position and velocity of a satellite. 


NUMERICAL RESULTS 


We now will give some results of an error analysis for a series of differ- 
ent cases for orbits about Mars and Venus , Moon, and Earth. However, let us first 
briefly describe the method of the error analysis. 

It has been shown that it is possible to write 

6m + D = A 6Y (0) (19) 


where 6m = a vector of angle deviations from the directions which would be 

computed if the primary and secondary were on the reference orbits, 

D = a vector which arises because the reference orbits were chosen to 
be Kepler orbits, 

A = a 2n x m matrix, n is the number of sightings of the secondary, and 
m is the number of unknowns in the problems^ 

6Y (0) = a 1 x m vector which represents deviations from their true values 
of the initial conditions and the multiplicative constants of the 
perturbing forces . 
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Equation (19) was derived by a linearization technique, and hence is only 
an approximation. However, we will assume errors in the measured angles to be 
10 seconds of arc. For these small errors, we may regard (19) as exact. 

Let us assume E(6m + D) = 0, i.e., the expected value of each of the 
measured angles is its true value. Also, the distributions of each measured 
angle are independent. And, finally, the standard deviation of each measure- 
ment is 10 seconds of arc. More precisely, for each i 

o = a( 6a^ cos e^) = a(6e^) = 10 seconds of arc. 

In general, 2n > m; i.e., there are more equations in (19) than unknowns. 
Since this is the case, a "least squares" solution of (19) may be sought. If 
the columns of A are linearly independent vectors, then the "least squares" 
solution of (19) yields 

E (6Y (0) ) = 0 

and 

E(6Y(0) 6Y(0)') = ct 2 (A'A)" 1 . 

Hence, the mean of the outputs will be their true values. However, even though 
the input errors are uncorrelated, the output error may be highly correlated. 

In any event, as a measure of the error in outputs due to the input errors, we 
will use the standard deviation of the output. 

The unit of length is chosen as the kilometer, and the unit of time is 
the hour. 


Mars 

For orbits about Mars, the secondary will first be chosen as one or both 
of its known natural satellites, Phobos and Deimos. Later, the effect of using 
an ejected probe will be studied. We assume Mars has an axis of dynamic symmetry, 
and hence we are justified in assuming a potential of the form (18) if the direc- 
tion of this axis is included in the unknowns. 
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Position and Velocity Only .- Let us first investigate the errors in position 
and velocity for the case in which the only force field is an inverse square 
central field; i.e., that is, no perturbing force exists. 

The orbital elements chosen for the spacecraft, Phobos, and Deimos are 
given in Table IV. 


Table IV 

Orbital Elements of Spacecraft, Phobos, and Deimos 


Body 

a (km) 

e 

i 

(JU 

0 

T (hrs) 
P 

Spacecraft 

4056 .4 

0 

45° 

0 

0 

0 

Phobos 

9400. 

0.017 

0.95° 

270° 

319° 

0 

Deimos 

23500. 

0.0028 

1.3° 

0 

305° 

0 


In Figure 6 is a plot of the distance between the spacecraft and Deimos, 
the spacecraft and Phobos, and the intervals of time over which these secon- 
daries will be obscured by Mars. 

Table V gives the resulting errors in position and velocity at t = 0. For 
this case, 132 sightings were uniformly spaced over 10 orbital periods of the 
satellite. Each sighting yields two equations, since both azimuth and elevation 
of the secondary were utilized. No sightings were utilized during the time inter- 
val in which the secondary was obscured by Mars. 

Table V gives the results for various cases. Given are the standard devia- 
tions of the position and velocity errors at t = 0. The components of these 
error vectors are resolved in the system R^(0); i.e., they are given in the 
radial, tangential, and cross-track directions (j = 1, 2). The 12 unknown case 
is that in which both the orbit of the spacecraft and Phobos are unknown. Two 
6 unknown cases are presented; one in which the orbit of Phobos is known, and 
sightings on it are used to compute only the orbit of the spacecraft. The 
second 6 unknown case utilizes Deimos instead of Phobos. Finally, an 18 un- 
known case is given in which sightings of both Phobos and Deimos are utilized 
to compute all the orbits of all three bodies. 
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CE BETWEEN VEHICLE AND SECONDARY (x|0° km) 


28 


DEIMOS 


1 6 i— 


K h 



PHOBOS 


OBSCURED INTERVALS 


ORBIT NUMBER 

4 5 6 7 8 9 /0 

II III II 


i i i i i i i » 

8 10 12 14 16 18 20 22 


HOURS 



Figure 6: Distance Between Spacecraft and Secondaries as a Function 

of Time 
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In all cases, the results indicate that the positions and velocities may 
be computed within a reasonable error. As might be expected, the case which 
yields the smallest errors is that in which an accurately known orbit of Phobos 
is utilized. Note that if the orbit of Phobos is treated as unknown, then 
the accuracy to which the spacecraft's position is determined is only slightly 
poorer . 

From this table, the overall conclusion is that the navigational problem 
is well conditioned. Moreover, the proper choice of problems is a 12 unknown 
problem in which the orbits of Phobos and the spacecraft are determined. 

Little is gained by utilizing sightings of Deimos. 

Table V was obtained by using a spacecraft orbital inclination of i = 45°. 
Several different inclinations were also utilized. The results were that the 
largest error was obtained for i = 0°, and the smallest for i = 90°. Hie 
spread, however, was not great. 


Table V 

Position and Velocity Errors at t = 0 
132 Sightings, 10 Orbital Periods 



Position Errors (km) 

Velocity Errors 
(km/hr) 

No. 

of 

Unknowns 

Body 

a( r ) 

a(s) 

a(z) 

2 

atf) 

a(5) 

a(z) 

s/c 

.044 

.141 

.081 

.168 

.380 

.126 

.246 

6 

Phobos known 

S/C 

.113 

.401 

.214 

.468 

1.06 

.330 

.583 

6 

Deimos known 

Phobos 

.071 

.190 

.079 

.217 

.116 

.059 

.060 


Deimos 

.375 

.432 

.167 

.596. 

.112 

.058 

.037 

18 

S/C 

.117 

.147 

.102 

.214 

.391 

.338 

.237 


Phobos 

.075 

.197 

.082 

.226 

.119 

.062 

.062 


S/C 

.125 

.159 

.112 

.231 

.420 

.362 

.262 

12 


2 = 's/d^ir) + a^(s) + (z) 
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r = error in radial direction 
s = error in tangent direction 
z = cross -track error 

Finally, in Figure 7 we plot what is defined as the total navigational 
error for the cases examined in Table V. We define total navigational error 
as , 

qT 

1 f V^O^t)) + <r 2 (s(t)) + a 2 (z(t)) dt 
(q-l)T 

9 . = 1 9 2 , «•» 

where T = orbital period of the spacecraft and all measurements are utilized 
such that 0 < t < qT. 

Unknown Potential Parameters and Coordinate Direction but Mass Known , -Let 
us now assume the mass of Mars is accurately known (in the next section, we 
allow for an unknown mass) and the potential is given by 

V(r, « - f=[- 1 +J 2 

We assume J^, J^, 5^, and are poorly known. The angles and g^ define 

the direction of the spin axis of Mars. 

In the previous section, we showed little was gained by sighting of Deimos 
(at least for the navigation problem and a low altitude spacecraft orbit). 

Thus, for the rest of our study of Mars, we will only utilize Phobos as the 
secondary body. The orbit of Phobos will be treated as unknown. The number 
of unknowns in the total problem is then 17 or 15 (in two cases we give results 
in which and g^ are known). 

Four different orbits of the spacecraft were investigated. These four 
orbits are defined in Table VI. For all orbits, (jd = Q = ^ =0. 
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Table VI 

Spacecraft Orbit 


Orbit 

No. 

a (km) 

e 

i 

Periapsis ht. 
(km) 

Apoapsis ht. 
(km) 

1 

4056.40 

0 

45° 

672.06 

672.06 

2 

4056.40 

0 

0 ° 

672.06 

672.06 

3 

3930.34 

0.114494 

0 

0 

100 

1000 . 

4 

8880.34 
1 

0.506737 

60° 

1000 

10000 . 


Values of the parameters were chosen as follows: 


§! = 317.9° 

S 2 = 35.3° 

J 2 = 2.011 x 10"3 

J, = J, =0 
3 4 , 

J 3 = -5 x 10 | 

J 4 = -4 x 10' 6 j 
No estimates of or J. are 


from [ 7 ] 
from [ 11 ] 

when orbits 1 and 2 were used 

when orbits 3 and 4 were used 
given in the literature. 


As expected, the position and velocity errors are now greater than those obtained 
for the case in which no perturbating forces exist (compare Table V with Table VII) . 
Also, for any particular case of the 17 unknown problem, a(z)>a(r)>a(s) . Deletion 
of the direction of the axis of symmetry as an unknown reduces a(z) and a(z) more 
than the other position and velocity error components, respectively. However, this 
deletion has little effect in reducing the potential parameter errors. 

At present, astronomic measurements yield £ 1 ( 6 ^) “10 a( 6 § = 0 ( 6 £ = 1°. Hence, 

the sighting can be expected to yield considerable improvement in our knowledge of 
the Martian potential. 

Unknown Mass. Potential Parameters, and Coordinate Direction . -Let us now con- 
sider exactly the same case as that presented by Orbit No. 3 in Table VII except 
we now assume that the mass of Mars is unknown. Hence, the problem now has 18 
unknowns. Results of this case are shown in Table VIII. 


From this table, we see the method is very poor if the mass of Mars and the 
orbit of Phobos are unknown. Note the large values of a(r) , a(s) , and a (6 m) . The 
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132 Sightings 



Position Errors at 
t = 0 
(km) 

Velocity Errors at 
t = 0 
(km/hr) 

Potential Parameter 
Errors x 10^ 

Axis 

Direction 
(secs, of arc) 


Orbit 

No. 

No. of 
Orbits 

Body 

a( r ) 

a(s) 

0(2) 

o(r) 

a(s) 

a(z) : 

ct( 6 J 2 ) 

o( 6 J 3 ) 


0 ( 6 ^) 

c (5f 2 ) 

No. of 
Unknowns 

1 

10 

Phobos 

S/C 

.08 

.14 

4.41 

.37 

19.54 

.24 

3.74 

1.03 

.06 

.43 

17.1 

35.4 

1.57 

.26 

7.15 

125 

622 

17 

1 

20 

Phobos 

S/C 

.07 

.14 

.81 

.13 

2.49 

.51 

.68 

.34 

.05 

.42 

2.28 

4.65 

1.49 

.12 

6.79 

25 

81 

17 

1 

10 

PhObos 

S/C 

.07 

.14 

.21 

.16 

.09 

.22 

.13 

.40 

.06 

.42 

.06 

.32 

1.57 

.23 

7.13 

--- 

— 

15 

1 

20 

Phobos 

S/C 

.06 

.14 

.21 

.12 

.08 

.19 

.12 

.30 

.05 

.42 

.07 

.30 

1.46 

.12 

6.65 

— 

— 

15 

2 ; 

10 

Phobos 

S/C 

.22 

.17 

.36 

.33 

5.61 
.17 ; 

.23 

.55 

.18 

.46 

4.1 

9.4 

3.66 

1.44 

35.7 j 

6 

165 

17 

2 i 

i 

20 

i 

Phobos 

S/C 

.11 

.15 

.27 

.29 

2.74 : 
.16 

.16 

.45 

.10 

.41 

2.0 i 

; 

1.40 

1.41 

8.81 

5 

81 

17 

! i 

. 3 ! 

j 

10 i 

Phobos 

S/C 

.10 

.16 

5.95 

.37 

29.4 
2 .90 1 

5.06 

1.35 

.07 

.54 

21.6 
55.9 1 

.57 

1.21 

1.05 I 

168 

863 

17 

1 1 

; 1 

'■ 3 ! 

20 

Phobos 

S/C 

.07 

.15 

1 .91 

; .28 

4 .49 
.53 

1 

.76 

.91 

.06 

.51 

3.37 

8.58 

.29 

.52 

.40 

28 

133 

17 

i 4 | 

10 

Phobos 

S/C 

.09 

; .10 

j .43 

I - 62 

5.59 

.40 

.32 

1.31 

.08 

.20 

4.0 

10.6 

j 

.65 

.36 

1.15 

16 

162 

17 

4 ! 

t 1 

;i i 

| 20 

Phobos 

S/C 

.08 

.09 

! .23 

j .44 

3.03 

.10 

.14 

.89 

; .06 
; • 18 

I 

*.17 

5.61 

i 

.24 

IL . . 

.14 

.55 

3 

88 

17 


Table VII: Position, Velocity, and Parameter Errors (Hass of Mars Known) 



Position Errors 
at t = 0 (km) 

Velocity Errors 
at t = 0 (km/hr) 

Potential Parameter 
Errors 

Axis 

Direction 
(sec. of arc) 

o 

• 

o 

55 

Body 

a(r) 

o(s) ct(z) 

a(r) ct(s) a(z) 


°(6J 3 ) 


a(6m) 

C7(6§ 1 ) 


Orbits 

m 

10 

Phobos 

S/C 

547 

206 

5.66 27.0 
.37 2.8 

4.8 464 20 

1.3 791 51 

1.97 x 10' 5 

1.26 x 10" 5 

3.33 x 10" 5 

.18 

160 

792 

20 

Phobos 

S/C 

358 

135 

.85 4.1 

.28 .50 

j 

.71 304 3.1 

.93 518 7.8 

1.23 x 10' 5 

5.61 x 10' 5 

2.06 x 10 -5 

.11 

25.8 

121 


Table VIII: Position, Velocity, and Parameter Errors (Mass of Mars Unknown) 









analytic reason for these poor results is obvious, for the equations of motion 
are of the form 


- G m r . - ,- N 
r = 5 - + g (r) 

Irl 3 


where g (r) is small if J^, an< * ^4 are sma H« 


Now let \R = r 
3 

\ M = m, where X is a constant. 

Since no lengths are measured, the factor X does not appear in the constraint 
equations if these equations are expressed in the new variables. However, in 
terms of the new variables, the equations of motion become 

r = g_m! +i 

PI 3 

Thus, we note, if g = 0, then the equations of motion are independent of \. 
Since g is small, it can be expected that r and m are nearly dependent, and 
thus cannot be well determined by the measurements. A penalty thus must be 
paid for not measuring lengths. 


In summary, sighting of Phobos from a spacecraft yields accurate determina 
tion of the spacecraft's orbit, the orbit of Phobos, and the potential para- 
meters. However, the mass of Mars cannot be treated as unknown in the total 
problem. 

Ejected Probe . -Thus far, we have utilized only Phobos as the secondary 
satellite in determining the potential of Mars. We have observed that if the 
orbit of Phobos and the mass of Mars are treated as unknown, then intolerably 
large errors are produced. This result occurs because no lengths are measured. 


Let us now investigate the total problem again utilizing Orbit No. 3 of 
Table VII for the case in which an ejected probe of known ejection speed is 
utilized. Note that the ejected speed introduces a known length. 
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Let the probe be ejected so that 
Av = 4 km/hr (ejection speed) 



The angles cy Q and e Q determine the direction of ejection as shown in Figure 8, 
We will assume these angles to be unknown, but a v known. 



Figure 8: Direction of Ejection with Respect to System 

which Moves with Spacecraft 

In Figures 9, 10, and 11 results of two cases are presented. One case 
utilizes 10 sightings per spacecraft orbital period on the ejected probe and 
10 sightings per period on Phobos; the other utilizes only the sightings on 
the probe. The first case in which Phobos is used contains six more unknowns 
than the second case. The two cases have 20 and 14 unknowns, respectively . 

Figure 9 gives the error in mass of Mars, navigational error of the 

spacecraft, and error in for both cases as a function of the number of 
orbital periods over which the data is gathered. Here, navigational error is 
defined as 
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ERRORS 



Figure 9: Total Navigation Error, Mass Error, and Error. 

Sighting of Probe Alone, and Probe and Phobos 
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Figure 10: Error in and ^ by Sighting of Probe Alone, 

and Probe and Phobos 
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ERRORS 



Figure 11: Errors in and Sighting of Probe Alone, 

and Probe and Phobos 
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y, = ( r ) + a 2 (s) + a 2 ( 2 ) 


at t = 0, 


and r, s, and z are defined in Figure 3. Note that the mass determination 
is always better if Phobos is used in conjunction with the probe. 


In Figure 10 we plot the error in and ^2 for t * ie tw ° cases# §2 t * ie 
codeclination of the axis of symmetry of Mars. Finally, in Figure 11 we plot 
the error in (right ascension) and J^. 

In general, there is no greatly significant difference between results 
obtained with and without sightings of Phobos along with the probe; except that 
in the first case, an accurate orbit of Phobos is determined. 

It is somewhat difficult to evaluate the improvement in accuracy afforded 
by sighting of the probe over that obtained by sighting Phobos alone from the 
data we have presented. Recall that Table VII (sighting of Phobos alone) was 
constructed by utilizing the assumption that the mass of Mars was known; how- 
ever, in Figures 9, 10, and 11 we assumed the mass to be unknown. Even though 
this is the case, the improvement obtained by the probe is about a factor of 
three for position, 15 for J^, 30 for J^, 23 for J^, and five for the direction 
of the axis of symmetry. 

If sightings of the probe and Phobos are utilized, the improvement in 
the orbit determination of Phobos (versus Phobos alone) is given in Table IX. 

Table IX 


Comparison of Accuracy of Position and Velocity at t = 0 of Phobos 
With and Without Probe Sightings, Orbit No. 3. Martian Mass Unknown With 

Probe, Known Without Probe. 


No. of 

No. of Measurements 

E 

WKGSSM 

t 


Orbits 

With 

Probe 

Without 

Probe 

With 

Probe 

Without 

Probe 

With 

Probe 

Without 

Probe 

10 

■SI 

264 

8.82 

1 

5.94 

22.2 

20 

MM 

264 

1.46 


0.84 

3.46 
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Venus 


^2 and Direction of Axis of Symmetry Unknown - Known Ejection Speed ."Of all 
the planets, Venus is the most nearly spherical. Its equatorial and polar radii 
are both estimated to be 6100 + 50 kilometers. Thus, it is difficult to deter- 
mine the coefficients of the higher order terms in the potential. No estimates 

are given in the literature. Here we assume all coefficients are zero except J^. 

The orbit chosen for the spacecraft is defined by 
a = 6587 kilometers 

e = .0151814 

i = 60° 

Q = (JL) = T =0. 

p 

These elements imply a periapsis and apoapsis altitude of 387 and 587 kilometers, 
respectively. 

We assume the spin axis to be coincident with the axis of dynamic symmetry 
of Venus. This assumption and data from [6] thus implies 

= 272.75° (right ascension of spin axis) 

^2 = 18.5 (co-declination). 

Venus has no known natural satellites, hence the secondary must be an 
ejected probe. We assume that the probe is ejected with known speed, but poorly 
known direction. It turns out that the total problem is unstable if the ejection 
speed is treated as unknown. We will discuss this point later. Also, we will 
show the effect of an error in the ejection .speed . 

We define the direction of ejection with respect to the coordinate R^(0) 
by two angles, <y Q and e Q , as pictured in Figure 8. 

Let us set 
Av = 4 km/hr 
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These parameters imply that the maximum distance between the spacecraft and 

-3 

probe in 10 orbital periods is 86.39 km, 86.50 km, and 86.51 km for J 9 = 10 , 

*-4-5 L 

10 , and 10 , respectively. Assume no errors in the mass of Venus and the 

ejection speed. There are 11 unknowns in the problem. These unknowns are given 

in Table X, along with results for various values of 

From Table X we note that as J 9 becomes smaller, the errors become larger. 

-6 ^ 

In fact, for < 10 the problem becomes so ill-conditioned that no solution 
could be found. The physical reason for this fact is that for small J^, the 
direction of the planet's axis of symmetry is poorly defined. It is undefined 
for J 2 = 0. 

For rapidly rotating planets (e.g. Jupiter and Saturn), can probably be 
estimated closely by assuming the planet to be in hydrostatic equilibrium. How- 
ever, for such a slowly rotating planet (sidereal period = 243.2 days retrograde), 
the effect of centrifugal force is second order. It is difficult for us to esti- 
mate for Venus, but it is highly likely that it is smaller than that of the 
Earth or moon. 

In Table XI we show the effect of adding the mass of Venus as an unknown. 
Recall that for the case of Mars, adding mass as an unknown greatly increased 
the errors. A significant, but not intolerable, increase in position and ve- 
locity errors are shown for the present case. This is because the ejection speed 
supplies a length measurement which was missing when only the direction of Phobos 
was measured. The accuracy of the mass determination, however, is not impressive. 
Our present knowledge (Table I) gives 6m/m = 1.5 x 10 

Note that (Table XI) deleting the direction of the axis of symmetry from 
the unknowns greatly decreases the errors in all remaining unknowns except 
(disappointingly) the mass. 

In Figure 12 we plot the errors as a function of the number of orbital 

-3 

periods over which data was gathered. Here = 10 , and 13 sightings/10 orbital 

periods were utilized. 


% 
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Table X: Errors in X. (0), J 2 , Direction of Axis of Symmetry, and Direction of Ejection: Speed of 

Ejection ana Mass of Venus Known 


Chosen 

Position Errors 

(km) 

Velocity 

Errors 

(km/hr) 

ct (6J 2 ) 

Axis Dir. 

(min) 

Ejection Dir. 

(min) 

J 2 

a(6x) 

a (6y) 

a(6z) 

G (6*) 

CT (6y) 

a(6z) 

x 10 6 

<7(6?]) 

o(6§ 2 ) 

cos e a( 6 sT) 

cr(6e) 

10' 3 

.042 

.451 

1.65 

4.63 

170 

98.1 

3.1 

.85 

27 

15 

17 

icf 4 

.041 

2.12 

6.60 

18.2 

1790 

1030 

3.3 

3.6 

280 

157 

183 

icf 5 

.041 

21.0 

62.9 

172 

19600 

11300 

3.6 

35 

3058 

1710 

2000 

10- 6 










— 


O 

T— 4 

— 

— 

-- 

-- 

-- 

-- 

-- 

-- 

-- 

-- 

— 


10 orbital periods of data - 130 sightings 


Table XI: Effect of Adding Mass as an Unknown. Also Effect of Deleting Direction of Axis of Symmetry 


Chosen 

J 2 

No. of 
Unknowns 

Posit ion 
Error 

2 (km) 

Speed 

Error 

£ (km/hr) 

g ( 6J 2 ) 
x 10 6 

Axis Direction 
(min) 

Ejection Direction 
(min) 

<?(6m) 

°(6? 1 ) 

a(6§ 2 ) 

cos e a(6a) 

a(6e) 

m 

IQ’ 3 

12 

30.9 

238.9 

3.2 

.90 

28 

17 

18 

.014 

10' 4 

12 

317 

2694 

3.4 

3 

294 

165 

192 

.185 

io' 3 

10 

29.4 

U5 

.77 

— 

— 

1 

1 

.013 

10- 4 

10 

37.8 

148 

.82 

— 

--- 

i 

1 i 

1 

.174 


10 orbital periods of data - 130 sightings 
2 = Va 2 ( 6x) + a 2 (6y) + a 2 ( 6z) 
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Figure 12: Effect of Additional Measurements and Orbits on 

Errors 
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In order to assure that the probe is not too far from the spacecraft, the 
angles of ejection for the case shown in Figure 12 were changed so that 

cy =0 
o 

e = 45 , where these angles are defined in Figure 8. 
o 

This ejection direction will be used in the rest of the analysis. Setting <* o = 0 
insures (at least to a first order) that no secular dependence is present in the 
distance between the probe and the spacecraft. During a period of 100 orbits 
the maximum distance between the probe and spacecraft is 3.91 kilometers. 

Note that even though more sightings were used for the case presented in 
Table XI than that presented in Figure 12, after 10 orbits the errors are smaller 
in the latter case. This results because a more judicious choice of the ejection 
direction yields a smaller distance between the spacecraft and probe. 

Several other cases were tried in which the planet* s mass was unknown, but 
errors in the mass determination were always quite large. At this point, we 
conclude that our method is somewhat poorly conditioned if the planet * s' mass is 
unknown In the rest of the analysis we assume the planet's mass is given. 

Let us look at one more case before leaving this section: 

(1) Mass known 

(2) a = 0 

o 

e = 45° 
o 

(3) Av = 4 km/hr 

(4) five sightings /orbit 

(5) J 2 = 10" 5 and l(f 6 . 

Conditions (1), (2), (3), and (4) will be used in the rest of our study of Venus. 
Results of this case are shown in Figures 13 and 14. Note that all errors are 
significantly larger for the smaller value of except <^(6^) itself. 

, ^ 21 ' anc * " 21 ^ n ^ nown -’-^ n order to overcome the singularity inherent in 
the above formulation for 3^=0, a different and more conventional formulation will 
now be made . 
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Assume a body has an axis of symmetry, and a potential which includes 
terms with Legendre polynomials up to degree k when the preferred coordinate 
system is utilized. Then as was observed in our discussion of the potential, 
a system slightly removed from the preferred system introduces first order 
Legendre functions if the removed system is utilized to express its potential. 


We now assume the potential of Venus to be 


V(r, <*, 0) 



So, C and are now used as unknowns instead of g ^ and Note that even 

if = 0, C 21 and S 21 have physical significance, but g^ and g^ do not. 


For the present study we set the direction of ejection so that a q = 0, 

= 45°. All orbital parameters are the same as before. However, five 
sightings are spread uniformly over an orbital period. Moreover, it is assumed 
the ejection speed is known as is 4 km/hr. The present problem has 11 unknowns. 

Results for the case J 2 - 10 C = 10 and = 0, are shown in 
Table XII. However, two other cases viz. = 10 C 21 = S 21 = 0; and 

= ^ = 0 yields almost identical results. This result is in sharp 

contrast with those obtained for the formulation in which the direction of the 
axis of symmetry was unknown instead of and Also, the error in initial 

velocity is much smaller in the present formulation. 


Table XII: Errors in Initial Position and Velocity, , ^21* anc * ^21 as a ^unc- 

tion of Number of Orbits, Mass known 


No. of 
Orbits 

£ (km) 

£ 

(km/hr) 

a (6 J 2 )xl0 6 

<J(6C 21 )x10 8 

c(6S 21 )x10 9 

5 

1.15 

1.37 

2.80 

28.7 

65.0 

10 

.82 

.95 

.98 

5.84 

11.3 

15 

.67 

.85 

.54 

2.19 

4.10 

20 

.58 

.67 

.35 

1.10 

2.03 

25 

.52 

.60 

.25 

.63 

1.14 

30 

.48 

.55 

.19 

.40 

.72 

35 

.44 

.51 

.15 

.27 

.49 

40 

.41 

.48 

.13 

.20 

.35 
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Effect of Ejection Speed Errors 


Extremely large errors in all parameters are obtained if the ejection 
speed is taken as an additional unknown in the total problem. The reason for 
this fact is that 

R 2 (t) - R 1 (t) = Av «S 2 (t) u q> (0) = 0, 

where = position of probe with respect to spacecraft, 

/\v = ejection speed, 

^2 - 3 x 3 matrix, and 

u = ejection direction, 
o 

Thus, the direction to the probe as a function of time is nearly independent of 
the ejection speed. Since only this direction is measured, the measurements 
weakly determine the ejection speed. 

The situation with respect to the direction of ejection is just the oppo- 
site. That is, the generic direction is strongly dependent on the ejection 
direction and thus the ejection direction is always accurately determined by 
the measurements. 

In all our previous analyses, we assumed the ejection speed was known 
exactly. More precisely, we assumed no systematic errors. The random errors 
in the measurement of the probe's direction then gives rise to random errors 
in the output parameters with mean zero and the reported standard deviation. 

Let us now investigate the effect of a systematic error in the ejection 
speed. In general, any systematic error will force upon us errors in the output 
parameters with non-zero means, but the standard deviation will be unaffected. 
Thus, our previous results may still be interpreted as the standard deviation of 
the parameter . 

Let us now give the systematic error produced by an error in the ejection 
speed for the five principal cases studied previously. (See Table XIII.) 

For these cases, we assume the mass of Venus known. 
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Figure 15 is a plot of the systematic error in the spacecraft's initial 
position and velocity, while Figure 16 yields the systematic errors in , 
and S 21 - These plots are for Cases III, IV, and V (Table XIII) and 
6Av * .01 Av (the error in the ejection speed is 1/100 of the true ejection 
speed) . We define the position and velocity errors plotted in Figure 15 
as the square root of the sum of the squares of the component errors . For 
the case shown the true magnitudes of the initial position and velocity are 
6,487 kilometers and 25,367 km/hr, respectively. 


From these figures, three main generalizations may be made as follows: 


(1) Smaller perturbations yield smaller systematic errors due to 
ejection speed errors. 

(2) In general (Case V in which the perturbing forces are zero 
offers exceptions) , the systematic errors in initial position, 
velocity, and C 21 increase with increasing number of measure- 
ments and measurement time interval. This is in contrast with 
the effect of random input error which always imply a decrease 
in output errors as more measurements and a longer time interval 
are used. 

(3) If the error in ejection speed is on the order of 1% or less, 
the resulting output errors are reasonably small. 


Table XIII: Cases for Which Systematic Error has Been Studied 


Case No. 

Parameters 

Random Error 
Results 

Comments 

i 

J 9 = io" 5 
§1> §2 

Figures 
13 and 14 

Systematic errors slightly larger 
in Case I than in II, but random 
errors were significantly larger 
in II than in I. 

II 

III 

j 2 = 10‘ 6 
§1> § 2 

j 2 = 10 -5 

C 21 = 10 ' 7 > S 21 = 0 

Table XII, 
Cases III, 
IV, and V 
yield almost 
identical 
results . 

Systematic errors smaller for 
smaller perturbations but random 
error nearly independent of magni- 
tude of perturbations. 

IV 

J 2 = 10 " 7 

C 21 = S 21 = 0 

V 

J 2 = C 2l 

J S 21= 0 
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SYSTEMATIC ERRORS 



Figure 15: Mean Error in Initial Position and Velocity as a 

Function of Number of Orbits. 6Av =.01 AA- 
Cases III, IV, and V 
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SYSTEMATIC ERRORS 



No OF ORBITS 


Figure 16: Mean Errors in J^, C^, an< * S 21 as a Function 

Number of Orbits. 6Av = .01 Av. Cases III, IV, 
and V 
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Finally, in Figures 17 and 18 we plot the value of |6A v |/Av which will 
produce a systematic error in any one parameter equal to the random error 
in that parameter. For example, from Figure 18 we see that for Case III 
after 20 orbits of measurements, the magnitude of the systematic error in 
the spacecrafts initial velocity is equal to the random error if |5Av| 

= .086 A v » If j 6 A v) < .086 Av, then the magnitude of the systematic error 
will be less than the random error. For Case IV, only the parameter C 
requires 1 5A v| < |a v | > for all other parameters 1 6Av[ > | A vj . That is, 
the error in the ejection speed is actually greater than this speed, to 
yield a systematic error in these other parameters equal to their respective 
random error. For Case V, no parameter requires j 6A v| < |A V | • 


Moon 

We assume the moon has a potential of the form 

v< r * *. 6 > --t[ 1+ I Z (f) 

j=2 i=0 

x cos i 9 + S_.^ sin i 0)J (20) 

where r, b , 0 = spherical coordinates of generic external point, 

G = universal gravitation constant, 
m = mass of the* moon, and 
r^ = equatorial radius = 1738.09 km. 

Also, we assume all forces on the spacecraft and probe are given by -VV. 

Three orbits about the moon will be utilized. The initial position 
and velocity and initial orbital elements are given in Table XIV. The 
maximum distance between the probe and spacecraft is dependent on the choice 
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Figure 17: Fractional Error in Ejection Speed Which Implies a Systematic Error 

in Output Parameters Equal in Magnitude to its Random Error. Case 
I and II. 
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Figure 18: Fractional Error in Ejection Speed which Implies a 

Systematic Error in Output Parameters Equal in 
Magnitude to its Random Error. Case III, IV, and V 




TABLE XIV: Initial Position, Velocity, and Orbital 

Elements of Orbits 1, 2, and 3 



Orbit 1 

Orbit 2 

Orbit 3 


s/c 

Probe 

s/c 

Probe 

s/c 

Probe 

x (km) 

2160 

2160 

1935 

1935 

1797.3 

1797.3 

y (km) 

0 

0 

0 

0 

0 

0 

z (km) 

0 

0 

0 

0 

0 

0 

x (km/hr) 

0 

2.828 

0 

2.828 

0 

2.828 

y (km/hr) 

2844.16 

2841.71 

3004.974 

3002.524 

3117.963 

3115.513 

z (km/hr) 

4926.23 

4927.65 

5204.767 

5206.182 

5400.470 

5401.884 

a (km) 

2400 

2400.00145 

— - 

2150 

2150.00116 

1997 

1997.001 

e 

.1 

. 100002 

.1 

. 100002 

.1 

.1 

i 

60° 

60.03° 

60° 

60.03° 

60° 

60.03° 

0 

0 

0 

0 

0 

0 

0 

0) 

0 

-.31° 

0 

-.29° 

0 

-.29° 

Tp (hr) 

0 

- .002 

0 

-.002 

0 

-.001 

Period (hr) 

2.930 

2.930 

2.484 

2.484 

2.224 

2.224 

hp (km) 

421.91 

421.91 

196.91 

196.91 

59.21 

59.21 

h fl (km) 

901.91 

901.92 

626.91 

626.92 

458.61 

458.61 


of C . . and S. . but over a time duration of 25 orbital periods, it is roughly 
Ji ji 

1.32 kilometers for Orbit 1 and 1.12 kilometers for Orbit 2. These distances 
are small because of a judicious choice of the ejection direction. In general, 
no secular dependence will be present in the distance between the probe and 
spacecraft if a zero component of velocity in the direction which is tangent 
to the spacecraft's orbital path is given to the probe. Thus, if a zero 
tangential component is achieved, the distance between the probe and spacecraft 
will be oscillatory and bounded. However, a non-zero component causes this 
distance to increase roughly linearly with time. For Orbits 1, 2, and 3 the 
probe was ejected so that equal velocity components were given in the radial 
and normal directions of the orbit and the ejection speed was 4 km/hr. 
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In Table XV, we present the four cases for which we will give results . 
The values of C^ and were taken from Reference [4]. In all cases we set 

C 21 " C 41 ~ C 42 " C 33 " C 43 " C 44 " S 32 " S 33 " S 44 * °* Ca8e 2 utilize8 the 
same orbit as Case 1. Cases 3 and 4 are the same as 2 except that I 3 utilizes 

a lower altitude orbit, while 4 utilizes the lowest altitude orbit (Table 

XIV). 


TABLE XV: Cases for Which Numerical Results will be Given 

(Multiply all values by 10"^) 


Constant 

Case 1 

Case 2 

Case 3 

Case 4 

C 20 

-2.07 

-2.07 

Same as 

Same as 

zu 

C 30 

.4461 

.4461 

Case 2 

Case 2 

C 40 

.2089 

.2089 



C 31 

.4346 

.4346 



C 22 

.2761 

.2761 



C 32 

0 

-.0522 



S 21 

-.4106 

-.4106 



S 31 

.1701 

0 



S 41 

-.1018 

-.1018 



S 42 

0 

-.0834 



S 43 

0 

-.0259 



Orbit No. 

1 

1 

2 

3 

Total No. of 

16 

18 

18 

18 

Unknowns 






The unknowns of the problem are the position and velocity of the space- 
craft at the time of ejection, the direction of ejection, and the potential 
parameters that are not identically zero in Table XV. 

For all results presented on the moon problem we chose to utilize 10 
sightings of the probe per orbital period of the spacecraft. Each sighting 
yields two equations. 
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Effect of Random Error in Direction 


Let us now investigate the effect of random error in Hi* measurement of 
the direction of the probe as seen from the spacecraft. These errors will 
imply an error in each of the outputs. 

Let us restate our major assumptions as follows: 


(1) No systematic error exists , and the only error source is random 
errors in each measured direction to the probe. 

(2) The expected value of each measured angle is its true value. 

(3) The error distributions of each measured angle are independent 
and identical with a standard deviation of 10 seconds of arc. 

(4) The random input errors are so small that each output error is a 
linear combination of the input errors. 

(5) Ten sightings per orbital period are made. Each sighting implies 
two equations. 


In Figure 19 we plot the initial position and velocity errors of the 
spacecraft due to the random input errors as a function of the number of 
orbits of the spacecraft. From this figure we note the following trends: 


(1) As the number of unknowns increases (Case 1 to Case 2) , the 
position and velocity errors increase. 

(2) As the orbital altitude is lowered (Case 2 to Case 3), the 
position errors decrease, but the velocity errors increase. 

In order not to crowd the graph, Case 4 (the lowest altitude 
orbit) is not shown, but its plot supports this conclusion. 

(3) The improvement in accuracy is rather gradual after 15 space- 
craft orbital periods. 


In Figures 20, 21, and 22 we plot the normalized errors in the parameters 
for Cases 1, 2, and 3, respectively. Here we define normalized error as 
Cj^ * ( <r C6 C j j ±1 , similar definitions for the S's. Results for 

Case 4 are not plotted, but a comparison of Cases 2, 3, and 4 is given in 
Table XVI. These three figures and the table yield the following generalization: 


(1) The potential parameter errors increase as the number of unknowns 
increases (Case 1 to Case 2) . 
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ERRORS 



5 10 15 20 25 


No OF ORBITS 


Figure 19: £ and £ as Function of Number of Orbits Over Which 

Data is Taken. Cases 1, 2 and 3 
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NORMALIZED ERROR IN FORCE PARAMETERS 



Figure 20: Normalized Error in Parameters as a Function of 

Data Interval. Case 1 
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TABLE XVI: Comparison of Results of Cases 2, 3, and 4. 

Data gathered over 25 Orbital Periods. 


Parameter 

Case 2 

Case 3 

Case 4 

£ (km) 

.142 

.131 

.124 

£ (km/hr) 

.121 

.131 

.137 

C 20 

.00270 

.00253 

.00250 

C 30 

.00633 

.00640 

.00630 

"40 

.147 

.0947 

.0727 

$31 

.117 

.0788 

.0567 

C 22 

.0491 

.0357 

.0299 

$21 

.0423 

.0285 

.0226 

C 32 

.0664 

.0515 

.0435 

*41 

.126 

.0725 

.0473 

^42 

.0113 

.00885 

.00680 

S 43 

.0250 

.0157 

.0118 


(2) As the orbital altitude is lowered (Case 2 to Case 3 to Case 4) , 
the parameter errors decrease. A significant decrease is obtained 
in all parameters except C^q and C^q (see Table XVI) . 

(3) Again the improvement is gradual .after 15 orbital periods. 

(4) The improvement in accuracy of is somewhat different from 

the other parameters in the interval from 5 to 15 orbital 
periods. Our conjecture for the cause of this more gradual 
improvement is the spacing of the sightings. Note that a factor 
of C- in (20) is cos 0. Now, the ejection was placed at 0 =0, 
and the probe nearly returns to the spacecraft for 0 = 2rrn, 

n = 1, 2, .... For this reason, no sightings were allowed within 
36° of 0 = 2rrn, which in turn generally implies a smaller coeffici- 
ent of than the other parameters at the time of probe sightings. 


In Table XVII, we compare a(6C^) and a(6S__) as obtained by our analysis 
(Case 3, 25 orbits) with those values reported by Lorell and Sjogren in 
Reference [ 12] . The values of the standard deviations from Reference [ 12] 
were estimated from the residuals obtained from actual data furnished by 
Lunar Orbiters I-IV. 
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TABLE XVII: 


Comparison of Standard Deviations of Parameters 
Obtained by Probe Sighting and Lunar Orbiter 
Analysis 

Standard Dev. x 10^ 


Parameter 

Probe Sighting 

Lunar Orbiter 

Probe/Orbiter 

C 20 

0.52 

1.43 

.36 

C 30 

0.29 

2.62 

.11 

C 40 

1.98 

1.90 

1.04 

C 31 

3.43 

0.25 

13.7 

C 22 

0.98 

2.49 

.40 

C 32 

.27 

0.58 

.46 

S 21 

1.17 

1.39 

.84 

S 41 

0.74 

0.51 

1.45 

S 42 

0.07 

0.35 

.20 

S 43 

0.04 

0.15 

.27 


From Table XVII, we see that our method compares fairly well with radar 
measurements. However, our error in is about 14 times larger. Moreover, 

we have not included the systematic errors which will arise in our method if 
an error in ejection speed is present. This point will now be considered. 

Effect of Error in Ejection Speed 

Thus far, we have assumed the ejection speed is given by Av = 4 km/hr 
and was error-free. Let us now assume the ejection speed is in error so that 
6Av = .04 km/hr. It seems reasonable that an ejection mechanism can be built 
so that | 6Av| < .01 Av. 

In general, if an error in ejection speed exists, a systematic error in 
each output will result. In Figures 23, 24, and 25 we plot the systematic 
error in the outputs as a function of time. Also shown are the "one-sigma" 
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PARAMETER ERROR x 10 









error bounds about the mean error. Case 3 is studied in these figures. 

From Figure 23 we note that the mean error in C^q increases as more data 
is utilized, but the variance about the mean due to the random error decreases. 
Similar remarks are true for the other parameters. In these figures we 
omitted the parameters C^, S 2i> anc * S 43* but t ^ ieir behavior is quite similar 
to that exhibited in Figure 25. 

EARTH ORBITS 

Thus far, the only perturbing forces we have considered were those which 
arise from a potential. With each harmonic in the potential, there was an 
associated unknown. 

In studying the application of our method for orbits about the earth, 
we will assume the potential is represented by three zonal harmonics such 
that the parameters , and are known beforehand. Additional perturb- 

ing forces due to drag and radiation pressure will be allowed. These non- 
conservative forces will introduce two unknowns. 

Two orbits will be used, a low altitude (632 kilometers) and a synchronous 
altitude orbit. The orbital parameters of the spacecraft are listed below. 

In all cases, five sightings per orbit will be utilized. 


Orbit Elements 

a 

e 

i 

0) 

0 

t 





p 

Low Altitude 

7,000 km 

0 

ON 

o 

o 

0° 

0° 

0 

High Altitude 

41,800 km 

0 

0° 

0° 

0° 

0 


Let us now briefly discuss the forces produced by drag and radiation 
pressure . 
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Drag 


The acceleration acting on a satellite caused by atmospheric drag may be 
approximated by 

a d = - (1/2 m) p v 2 AjjCp v 

density of the atmosphere at the satellite, 
cross-sectional area of satellite perpendicular to v, 
drag coefficient, 
air speed, 

unit vector in direction of satellite velocity with respect to atmos- 
phere , and 

mass of satellite. 

If we assume the atmosphere rotates with the Earth as a rigid body, then 
v v = R - | R| Q cos 6 e, 

where Q = angular rotation rate of the Earth about its axis = 7.29215506 x 10 ^ 
rad/sec , 

5 = declination of satellite, 
e = unit vector in local East direction, 

R = position vector of satellite. 

The drag coefficient, C^, is not constant, but is a function of angle of 
attack, temperature, and velocity. A brief, but enlightening, discussion of 
this parameter is given in Reference [13] . 

Here we will assume ” P is constant, but unknown. Hence, 

a d = - (<= d Ap v 2 /2 m) v. 


where p = 
^ = 



v = 

V = 


Radiation Pressure 

The acceleration of a satellite due to radiation pressure is given by 

a = - 1/m v v P A s, 

p Y p * 
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where y = a factor which depends on the reflecting properties of the satellite, 
= 1 + ,5k 9 for a diffuse reflector, 

= 2-2\/3, for a specular reflector; here 
\ is the albedo of the satellite, 
v = 1, if satellite is sunlit, 

= 0, if not sunlit, 

P = solar radiation pressure in the vicinity of the Earth (approximately 
4.5 x 10’ dynes/cni ), 

Ap = cross-sectional area of satellite perpendicular to s, and 
s = unit vector in direction of sun. 

Here we will assume P is constant, but unknown. Nominal parameters for drag 
and radiation pressure chosen here are listed in Table XVIII. 


TABLE XVIII: Drag and Radiation Pressure Parameters 


Parameter 

Spacecraft 

Probe 

m 

3000 kg 

.3 kg 

P 

4.5 x 10 (dynes /cm 2 ) 

4.5 x 10 5 (dynes /cm 2 ) 


„ A 2 

, 2 

II 

8 x 10 cm 

40 cm 

C D 

2.2 

-2.2 

Reflector 

specular 

specular 

\ (albedo) 

0.5 

0.9 

Y p A n / m 

10 v cm/sec 

— 6 2 

4.2 x 10 v cm/sec 

P 

p (low alt.) 

5.5 x 10 *kg/cm^ 

5.5 x 10" 16 g/cm 3 

p (high alt.) 

0 

0 

P ^ 
2m (10W) 

1.613 x 10 _17 /cm 

1.473 x 10* 16 /cm 
! 


The values of mass and area in Table XVIII imply that if the spacecraft and 
probe were spherical, then 


spacecraft radius 
probe radius 
mean spacecraft density 
mean probe density 


1. 12 meters , 

2. 52 cm = 1 inch 

3 

.498 g/cm , and 
.141 g/cm^. 
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Numerical Results for Earth Orbits 


For our first case, let us consider the low altitude (a = 7000 kilometers) 
orbit. As with the case of the Moon, let the probe be ejected with no tan- 
gential component of velocity. Also, let £v = 4 km/hr and 5 /^v = .04 km/hr. 
This error in ejection speed will produce a systematic error in each output 
parameter. Moreover, a random error will also exist because of random errors 
in the measurement of the probe 1 s direction. As before, these random errors in 
direction will be chosen with a standard deviation of 10 seconds of arc. 

Our present problem has 10 unknowns; six define the initial position and 
velocity of the spacecraft, two define the radiation and drag parameters, and 
finally, two are parasitic and define the probe ejection direction. 

Five sightings per orbital period of the spacecraft will be used. Each 
sighting will yield two equations. 

Figure 26 is a plot of the mean error in initial position and the mean 
error in the initial velocity as a function of the number of orbits over which 
the measurements are taken. Also shown are the one sigma error bounds about 
the mean. These errors about the mean are produced by random input errors. 

Note that the random errors are almost insignificant compared to the systematic 
errors. Also, the systematic error in initial velocity increases quite 
rapidly as data over a longer interval is utilized. 

In Figure 27 we plot the normalized systematic error in the radiation 

pressure and drag parameters. The random error in these parameters is about 

-2 

10 times the systematic error and hence too small to be shown in this figure. 

Note the error in the radiation pressure parameter, 6P, is almost as large as 

-3 

P. However, for the drag parameter 6c^ = 10 c^ so the measurements will 
yield a good estimate of the product p C^. Note that only the product is 
calculable and not the individual factors. 

It is not surprising that the 6P/P » 6c^/c^, for the acceleration due to 
drag is about 20 times greater than that due to radiation pressure. Moreover, 
no effect is produced by solar radiation over about 1/3 of each orbit. 
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POSITION AND VELOCITY ERRORS 






Since the effect of a systematic error in ejection speed outweighs the 
effect of random measurement errors in the previous case, it is highly de- 
sirable to reduce or eliminate the ejection speed error. We propose to ac- 
complish this elimination by releasing the probe with zero differential speed, 
i.e., a v = 6Av = 0. Drag and radiation pressure then is allowed to separate 
the spacecraft and probe. It is argumentative as to whether or not such a 
release is possible, but we will assume it is here. 

The total problem has now only eight unknowns, for we lose the two un- 
knowns associated with the ejection direction. 

In Figure 28 we then plot the errors in initial position and velocity due 
to the random errors. Two orbital cases are shown: the low altitude case 
and the high altitude (synchronous). Note that the errors decrease with more 
measurements and are considerably less than those shown in Figure 26. Also 
note the plot begins at 10 orbits instead of five orbits. This is because no 
measurement was allowed if the distance between the probe and spacecraft was 
less than 100 meters. If this restriction were not imposed, parallax may 
produce an additional error. After five orbits, the distance between the probe 
and spacecraft for the low and high orbits is 420 and 910 meters, respectively. 
In Figure 29, we plot the normalized error in radiation and drag parameter for 
the low altitude orbit, and the noramlized error in the radiation parameter for 
the high altitude orbit. It was assumed c^ = 0 for the high altitude orbit. 

Note that for the low altitude orbit the normalized drag error is about 

-2 

10 times the normalized radiation pressure error. But, the normalized 
radiation pressure error for the high altitude orbit is smaller than either 
error in the low altitude case. 

In conclusion, it appears that probe sightings are quite effective in 
calculating the parameters associated with drag and radiation pressure if 
systematic error in ejection speed can be avoided. 

As a final note of interest, we plot the distance between the probe and 
spacecraft as a function of time for the three cases. We offer the following 
explanation for the somewhat strange behavior shown in Figure 30. 



POSITION AND VELOCITY ERRORS 



10 15 20 25 30 35 40 


No OF ORBITS 


Figure 28: Initial Position and Velocity Errors for Low and 

High Altitude Earth Orbits. 6Av = Av = 0 
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Figure 29: Normalized Error in Drag and Radiation Parameters 

Caused by Random Measurement Errors. 6Av = Av = 



DISTANCE BETWEEN PROBE 8 SPACECRAFT ( km ) 


30 



0 4 8 12 16 20 24 28 32 36 40 

No OF ORBITS 

Distance Between Probe and Spacecraft as a Function of Time. 


Figure 30: 


Three Cases 



(1) For the high altitude orbit in which the ejection speed is zero 
(Av *= 0) , the only effect which separates the probe and space- 
craft is the differential radiation pressure. But radiation 
pressure is approximately in the direction of motion over 
one-half the orbit and opposite the direct over the other 
one-half orbit. Hence, the distance may be expected to be 
oscillatory. 

(2) For the low altitude orbit in which Av, the main effect which 
causes separation is drag. Drag always acts to oppose the 
motion. Hence, the separation distance monotonically increases. 

(3) If no drag or radiation forces exist, Av ^ 0 but small, and 
ejection is such that no orbital tangential component of 
velocity is given to the particle, then the separation distance 
would be oscillatory. If drag and radiation are added, the 
effect shown in Figure 30 for Av «= 4 km/hr is reasonable. 
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CONCLUDING REMARKS 


The major results of computer simulations indicating the manner in which 
instrument measurement errors propagate through the system into errors in the 
computed quantities is given below. 

Mars Navigation Problem (Natural Satellites) 

1. For the assumed conditions navigational accuracies ranging from 0.2 to 
1.0 kilometer are obtained when viewing either Phobos or Deimos or both. 

2. Introducing unknowns associated with the gravitational potential does 
not significantly degrade the accuracy of the spacecraft orbit determination; 
however, it does degrade the orbit determination of the artificial satellites 
(compared with the case in which the perturbing parameters are known). 

3. Approximately two or three orbits of observations are needed to get the 
navigational error down to the region where the improvement is gradual and 
roughly proportional to the reciprocal of the square root of the number of ob- 
servations. (Past studies have shown that one orbit is sufficient to arrive at 
this condition when using a test probe system. The fact that more orbits are 
taken here is a reflection of the fact that the orbital periods of the bodies 
being viewed are several times longer than that of the spacecraft orbital period.) 

4. As a navigational system, the accuracy when viewing Deimos is poorer 
than when viewing Phobos . One of the reasons for this is that Deimos is simply 
further away and the angle errors propagate in proportion to the length of the 
sight line. Another reason is that Deimos 1 orbital period is longer and for a 
given observation time its motion around the central body is smaller. 

Mars Perturbing Parameter Problem (Natural Satellites) 

J 2 3 t * ie number of unknowns is expanded to include three terms 

in the gravitational potential (making a total of 15), a well conditioned prob- 
- 7 - 4 

lem exists: 10 < a( 6X) <10 , i = 2, 3, and 4. 

* Only the unknown parameters are listed. 
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J 2 , J^, , g^, g 2 . -Adding the unknown direction of the spin axis, and 

g 2 , creates a problem in 17 unknowns. This problem is also well conditioned and 
the accuracy of the computed direction of the axis of dynamic symmetry varies 
from 1.0 to 10 minutes of arc. Some of the properties of this type, of system are; 

a. The terms in the gravitational potential only slowly improve in 
accuracy as the number of spacecraft orbits over which measurements 
are made is increased from 10 to 20; 

b. The spin axis direction determination is improved by a factor of two 
to eight by extending the sampling period from 10 to 20 orbits; 

c. An inclination of 60° generally produces a solution which is a factor 
of 10 more accurate than that produced by an inclination of 0; and 

d. Changing the spacecraft altitude is not a dominant factor in changing 
the accuracy with which the unknowns can be computed. 

J 2 , J^, J^, g^, g 2 , m. -Adding the planetary mass as unknown expands the 
problem to a total of 18 unknowns. This produces a system which: 

a. Errors in the computed spacecraft position become very large; 

b. The error in the mass determination is also large; 

c. Interestingly, there is no significant change in the accuracy with which 
the other parameters are computed (compared with the case in which mass 
was assumed to be known). 

d. The addition of a linear measurement (such as would be supplied by radar) 
will rectify the above deficiency of systems in which mass is an un- 
known and permit an accurate computation of mass to be made. 

Mars Ejected Probe (J 2 , J 3 , , g^ g 2 , m) 

Probe and Phobos . -Introducing an ejected probe expands the problem to a total 
of 20 unknowns in that two unknowns associated with the ejection direction must be 
added. The ejection speed cannot be used as an unknown without ill-conditioning 
the problem. It is assumed this speed is known exactly. 

A significant improvement is obtained by use of the probe in all parameters. 
Moreover, the orbit of Phobos may be found more accurately by about a factor of 
four. 
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Probe Alone . -Deleting sightings of Phobos reduces the probelm to 14 unknowns. 
No great change in the accuracy of the parameters is obtained except that the 
error in m is about a factor of two larger. 

Venus (6 A V s 0) 

J 2* 5]_» §2' -As J 2 becomes smaller, the errors become larger. The problem 
is ill-conditioned for = 0. 

J 2 , § 2 > m. -Adding the mass of Venus as an unknown does not significantly 

degrade the determination of and §£> but does degrade the initial posi- 

tion and velocity. The accuracy of the mass determination is poor and becomes 
poorer as becomes smaller. Mass should not be inserted into the problem as 
an unknown. 

J 2’ C 21 ’ S 21'” In contrast to case 1 C J 2» §2^ in wbicb §i and §2 were un “ 

known instead of C 2 ^ and S 2 ^ , the errors are almost independent of the magnitude 
of C 2 i , and S 2i * Tbe sin S u l arit y is thus removed by this formulation. 

Venus (6Av ^ 0) 

If an error exists in the ejection speed, a systematic error in each output 
will result. In general, as the perturbing forces become smaller, the systematic 
errors become smaller. However, as the time interval over which the measure- 
ments are taken is increased, the systematic errors are increased. But, in any 
event, an error in the ejection speed produced a very small error in the outputs. 

J 2 9 ?i» § 2 ' " Consi der the ratio of the systematic error in a quantity (pro- 
duced by an error in ejection speed) to the random error in that quantity 
(produced by random errors in measuring the direction to the probe). This ratio 
is relatively large for the initial position and § ^ and small for other outputs. 

If |6Av| < .01 Av, then the systematic error in all outputs is relatively 

small if measurements are made over less than 12 orbits and J = 10 . If 

- 6 ^ 

J 2 = 10 , then they are relatively small for 32 orbits. 
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J 2* C 2l > S 21*” For t ^ lis set unknowns, C ^ has the largest relative error, 
but if < 10“^ the systematic errors are relatively small. 


Moon 


1. The random errors in initial position and velocity are less than 0.2 kilo- 
meter and 0.2 km/hr, respectively, after 15 orbits of data collection 

2. A higher altitude orbit yields larger errors in the harmonics, except 
for C 2Q and C 3Q . 

3. The errors in the harmonics due to the random input errors have normalized 
values roughly between 0.1 and 1 after five orbits of data collection, 

and between 0.01 and 0.1 after 20 orbits. 

4. If an error is present in the ejection speed, a systematic error in each 
harmonic will result whose properties are as follows: 

(a) As a longer data collection time interval is utilized, a larger 
systematic error generally occurs. 

(b) For some of the harmonics, the systematic error is about 10 to 20 
times larger than the random error if 6Av = .01 Av. 


Earth 


1 . 


2 . 

3. 


Depending on the error in ejection speed, the product of atmospheric density 
(at the satellite) and the drag coefficient, pC^, can be determined to 
between one part in 1,000 and 10,000 by utilizing a low altitude orbit and 
a two inch diameter probe. 

If there exists an error in ejection speed of <5Av * .01 Av, then the error 


in the radiation pressure, 6P, is nearly as large as P. 

If the probe is released so that 6Av = Av ■ 0, then the following is true: 

(a) For a low altitude orbit, the drag forces are much greater than the 
solar radiation forces, and, hence, the normalized error in the 
radiation pressure is about 100 times greater than that of the drag 
parameter . 

(b) For a high altitude orbit, the drag forces are negligible. After 30 
orbits of data, the normalized radiation pressure error, 6 P/P, is 
about 2 x 10 -5 . 
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APPENDIX A 

ACCURACY OF ORBIT DETERMINATION 

We have proposed and utilized a method of orbit determination which relies 
on computing position and velocity by a linearization with respect to a reference 
trajectory. The reference trajectory is chosen to be a set of discontinuous 
Kepler arcs. 

In order to test the accuracy of our method, the sponsor supplied a computer 
program called Lungfish. This program computes the position by a numerical 
method of solution to the differential equation. The three second order differ- 
ential equations are used to obtain a 12th order Runge-Kutta type expansion for 
the solution. 

In Figure 31, we plot the difference between the Lungfish solution and our 
solution as a function of time. The orbit in question is a low altitude orbit 
about Mars (Orbit No. 3, Table VI). Only two components are shown, but the 
third component is similar to those shown. We feel that curves shown can be 
interpreted as more simply the difference between the two solutions and is the 
error in our solution. We conclude this because the Lungfish solution implied 
a nearly constant kinetic energy and component of angular momentum which is 
parallel to the symmetric axis (within 10 ^) . 

Note the maximum error in our solution is 270 meters after 20 orbits. 

It is felt that this error has a negligible effect on the results of our error 

analysis. A test of the accuracy of our solution was also made using an orbit 

-4 

about Venus with = 10 . In this case, the maximum error was only five 

meters. It is not surprising that the error in our solution is considerably 
smaller for smaller perturbations; for our solution is exact for zero perturba- 
tions . 


APPENDIX B 

USE OF PROBE FOR AN ORBIT ABOUT SATURN 
In our study of orbits about Venus, we found that the measurement of an 
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Figure 31; Error in X and Z - Components of Position Due to Method of Solution. Orbit about Mars. 



ejected probe's direction as a function of time yielded a poor determination of 
mass and a relatively poor determination of the direction of the axis of symmetry. 
We reasoned that these facts were caused because of the almost spherical 
character of Venus. To test our reasoning, we now study an orbit about Saturn 
which is a highly oblate planet. 

*4 

We now assume J = .0148, = 0, and = -3.5 x 10 , and an orbit whose 

semi-major axis has length 20 % larger than the planet's equatorial radius, and 
eccentricity = .01. The following results are obtained for no error in 
ejection speed. 

TABLE XIX: Standard Deviation of Parameters for an Orbit About Saturn. 

Mass Unknown 



< 4-4 

o 

o 

a 

Orbits 

Parameter 

10 

20 


-4 

-4 

a (6 m) /m 

3.9 x 10 

1.9x10 

CT (6? 1 )(sec of arc) 

43 

17 

<j(6g 0 )(sec of arc) 

20 

6 

z 

-6 

-7 

a(6 J 2 ) 

2.4 x 10 

-6 

7.3 x 10 

-7 

o(6 J,) 

1.5 x 10 ° 

8.0 x 10 


-6 

-7 

ct ( 6 j 4 ) 

1.0 x 10 

3.3 x 10 


Except for errors in , the errors shown in Table XIX are somewhat less than 
those shown in Figure 12. However, the determination of the mass is still not 
impressive . 

APPENDIX C 

CHANGE IN DIRECTION PRODUCED BY HARMONICS IN MOON'S POTENTIAL 

It is of interest to plot the change in direction of the probe produced 
by the small perturbing forces. We will do so now, utilizing a Moon orbit 
given by Orbit 1, Table XIV. To compute this change in direction, initial 
conditions and a Kepler field were first chosen. The right ascension of the probe 
was then computed as a function of time. The same initial conditions were 
rechosen but perturbing forces were added to the Kepler force. A new right 
ascension as a function of time was computed. The difference of these right 
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ascensions is shown in figure 32 (perturbed- -unperturbed) . Two cases 
are considered: one in which the oblateness term C^q produces the 

perturbing force, and a second in which this force is produced by the 
harmonics of Case 1, Table XV* 

For the cases shown, the nominal orbital period is 2*93 hours. 
Also, the ejection direction was chosen so that the probe nearly re- 
turns after each orbital period. Because a small change in position 
produces a large direction change when the probe is near the space- 
craft, large («* + 4°) direction changes are apparent near multiples 
of the orbital period. Also, the deviation in azimuth as pictured 
in figure 32 is large enough so that our intuition is not strained by 
the fact that 10 arc second errors are small enough to yield results 
of useful accuracy. 

The deviation in elevation is not shown, but this angle deviation 
is similar to that shown in the figure. 
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DEVIATION IN AZIMUTH PRODUCED BY HARMONICS 
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